Magnetic Brane-worlds 



John D. Barrow* and Sigbj0rn Hervik^ 

DAMTP, 
Centre for Mathematical Sciences, 
Cambridge University 
Wilberforce Rd. 
Cambridge CB3 OWA, UK 

February 7, 2008 



Abstract 

We investigate brane-worlds with a pure magnetic field and a perfect 
fluid. We extend earlier work to brane-worlds, and find new properties of 
the Bianchi type I brane-world. We find new asymptotic behaviours on 
approach to the singularity and classify the critical points of the dynamical 
phase space. It is known that the Einstein equations for the magnetic 
Bianchi type I models are in general oscillatory and are believed to be 
chaotic, but in the brane-world model this chaotic behaviour does not 
seem to be possible. 



1 Introduction 

Developments in string theory have inspired the construction of brane- worlds, 
in which standard-model gauge fields are confined to our three-dimensional 
brane-world, while gravity propagates in all the spatial dimensions. A sim- 
ple five-dimensional class of such models allows for a non-compact extra di- 
mension via a novel mechanism for localization of gravity around the brane at 
low energies Q ||. This mechanism is the warping of the metric by a negative 
five-dimensional cosmological constant. These models have been generalized to 
admit cosmological branes ||, and they provide an interesting arena in which 
to impose cosmological tests on extra-dimensional generalizations of Einstein's 
theory. One of the unusual features of this structure is that brane-worlds can feel 
the effects of non-local anisotropic stresses in the bulk dimensions. These gravi- 
ton stresses are unconstrained by dynamical equations specified on the brane but 
for a wide range of behaviours they completely determine the evolution of small 
anisotropics on the brane-world and can slow the their decay to such an extent 
that their effects are directly detectable today in the temperature anisotropy 
of the microwave background radiation^. The evolution of anisotropic brane- 
worlds has been studied recently by several authors ]|, ^, fj], ||, ||, They 
consider only the behaviour of simple anisotropics in universes with perfect fluid 
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matter sources. However, it is known that the evolution of simple anisotropic 
universes with isotropic three-curvature is very sensitive to the presence of trace- 
free matter sources with anisotropic pressures, O, |l3). Such stresses are 
inevitable in anisotropic cosmologies at early times because of the presence of 
collisionless gravitons, collisionlcss asymptotically-free particles, and electric or 
magnetic fields. These stresses also mimic the effects of more general anisotropic 
curvature anisotropics in more general anisotropic universe (because the latter 
can be viewed as an effective stress of long- wavelength gravitational waves). 
Small anisotropic pressures serve to slow the decay of the shear anisotropy to a 
logarithm of the time during the radiation era and to a slow power law during 
the dust era. These anisotropic stresses can arise because of intrinsic anisotropic 
stresses on the brane or via the induced graviton stresses from the bulk. In this 
paper we are going to consider the detailed behaviour of the simplest example: 
a pure magnetic field on a flat anisotropic brane with Bianchi type I geometry. 
Whereas the study of ref. Q was interested in the late-time behaviour we shall 
be interested primarily in the behaviour as t — > 0. 

The possibility of isotropization and homogenisation of cosmological models 
have been an important study since the original equations for gravitational inter- 
action where put down in 1915 by Einstein. The general relativistic anisotropic 
models called the Bianchi models have been studied in great detail both in gen- 
eral relativistic cosmologies [[l4| O, [ll| and to some extent in the brane-world 
scenario . A feature that has been known to exist in the Bianchi models are 
oscillatory solutions JO}, 0, and even chaotic solutions fl|, EJ|, |(| |2l], |j| || . 
The chaotic behaviour in the general relativistic vacuum Bianchi type VIII and 
IX models has been well studied, and chaotic behaviour for the general rela- 
tivistic Bianchi type I model with a magnetic field has been conjectured jl6| 
and studied for the case of form fields . There are many similarities between 
magnetic Bianchi universes and the behaviour of pre-big-bang string cosmolo- 
gies |25| . For further discussion of magnetic field effects in general cosmological 
models see also refs. pBL p7| . 

In the brane-worlds model the dynamics change because of the influence of 
terms varying as the square of the isotropic and anisotropic pressure as t —* 0. 
We investigate the Bianchi type I brane-world with a perfect fluid and a pure 
magnetic field in the context of a possible chaotic behaviour. This extends 
the general relativistic study by LeBlanc jl6| to the brane-world scenario. We 
will show that unlike the general relativity case, chaotic behaviour is probably 
completely absent for all perfect fluids as a result of the non-linear stresses. 

This paper is organized as follows: First we put down the fundamental equa- 
tions of motion. We use the representation of LeBlanc to facilitate comparisons, 
and introduce expansion- normalized variables. Since we are primarily interested 
in the initial singularity, we find some approximate solutions near the initial sin- 
gularity for different magnetic configurations. All of these approximate solutions 
correspond to equilibrium points for the full set of equations, and so we provide 
a stability analysis of the equilibrium points found. 
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2 Equations of motion for a Bianchi type I brane- 
world 



The Bianchi cosmological models can be classified according to the commutation 
relations of their invariant 1-forms, aA We will use an orthonormal frame of 
vector fields {e^}, and align the frame so that eo = d t , the fluid vector field. 
This implies that the spatial triad {e a } spans the tangent space at each point 
of the group orbits. Using the commutation relations for this vector basis, the 
commutator functions C a ^(t) are given by 

[e a ,ep] = C 7 a/ 3(t)e 7 (1) 

In the Bianchi Class A case we have the relations^] 

C a ob (t) = -® a b (t) - e a bc Q c (t), C\ c = e bcd n ad (2) 

The Q%{t) are the frame components of the volume expansion tensor of the 
fluid, or the extrinsic curvature. The vector f2 c (i) is the angular velocity of 
the spatial triad with respect to a set of Fermi-propagated axes. The extrinsic 
curvature is separated into its trace, 0, the volume expansion, and its 
trace-free part a M „, the shear: 



®nv = -^Qh^ + (3) 



We will also assume that we have a non-zero electromagnetic (EM) field and a 
non-interacting perfect fluid. The perfect-fluid energy momentum tensor is 

T lv = (p + P)u^u v + pg^u (4) 

Further, we will assume a barotropic equation of state for the fluid for the 
pressure p and density p: 

p=(j-l)p (5) 

where 7 € [0, 2] is constant. The perfect fluid solves the usual conservation 
equation 

p + 76/9 = (6) 

The energy-momentum tensor of the EM field can be expressed as 

T ^ M = PEMU^v+PEMh^ + n^, (7) 

where Pem an d Pem are the energy density and the isotropic pressure, re- 
spectively, and the tensor h^u — + u^Uu projects orthogonal to the on 
the brane, and is the tracefree anisotropic pressure tensor. In the case of 
an electric field E and magnetic field B the energy-momentum tensor can be 
written in general as 

T ^ M — --^F^Fapg^ + F^ a F v a (8) 



1 Latin indices have the range 1-3, while Greek indices have the range 0-3. 
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where the anisotropic stress is 



-E I1 E„-B I1 B V + -{E 2 + B 2 )h^ v (9) 



and the pressure and density are given by 

p EM =?>PEM = \(E 2 + B 2 ) (10) 

Demanding that the perfect fluid and the EM field be separately conserved 
requires eqn. (||) 

Pem + ®(Pem + Pem) + a^n^ = 0. (11) 

For a Bianchi type I model the commutator between the vectors of the spatial 
triad all commute so 

C a bc = (12) 

This allows us to have a pure magnetic field |28|. This is not the case for all 
Bianchi models. We shall therefore specialise to the case of a Bianchi type I 
model with a vanishing electric field. 

We now consider the evolution equations on the brane |2j], for the expan- 
sion scale factor, shear, and energy densities. The latter consist of the perfect 
fluid density, the magnetic energy and the dark energy term U. The non-local 
energy conservation equation for the dark energy U is 



u + -eu + a^V^ 



a: 4 



Y2 ^ V + HP + PK"*> + evr^vr^ - a^ir atl wj*] (13) 

where the overdot means u^V^ = d/dt, is the unconstrained bulk graviton 
stress, and k 2 = 8irG. The constraint equation is 



H 2 = ^ + ^^-^TZ 
3 6 6 



pi 



2\J 4A 



2U 



, , ( 14 ) 



where 7Z is the 3-curvature of the brane, A is the brane tension, A the cos- 
mological constant, and H = ^0 is the Hubble expansion parameter. The 
Gauss-Codazzi shear propagation equations are: 



6 

<?<nv> + = ^^^v + 77T [-(ft + ^p)^ixu + ir a<fl Tr a u> ] + —^tV^u - TZ <fl u> (15) 

Z A K A 



where A <uu> for a tensor A^u, means the projected, symmetric and trace- free 
part [ p9| |30[ | . The volume expansion propagation is governed by Raychaudhuri's 
equation 

6 + le 2 + a^a^ + X -K 2 (p + 3p) - A = 

' [n i (2p 2 +3pp) + 12U] (16) 



2\n 2 
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Here, for simplicity we shall assume that Vp. v = in order for the equations to 
close and we set 7?.< M „> = in the Bianchi I universe and also assume A = 0. 
The dark energy, 11 (which can in principle be of either sign) on the other hand 
cannot be set to zero in the presence of a non-zero anisotropic stress. In addition 
to these equations the magnetic field satisfies Maxwell's equations: 

B ll = ~QB ll + a llv B u (17) 

Since we shall be concerned with the early time behaviour we can ask if there 
is an initial singularity. In the brane-world we can write the energy-momentum 
tensor as a sum of the classical energy-momentum tensor and a term that come 
from brane-effects only. A sufficient condition for a singularity is that the clas- 
sical energy-momentum tensor obeys the strong energy condition (SEC), and 
that U > 0. We must however emphasize that we can still have a singularity 
even if this condition is violated. 

The equations in the orthonormal frame reduce to the following first-order 
system: 

e = a t e (is) 

B a = d t B a -e abc n c B b (19) 
&ab = d t a ab + 2a^ a e b ^ cd fl c (20) 
n a b = d t n ab + 2TT^ a e b ) cd Vl c (21) 

We still have an unused gauge freedom to choose the orientation of the spatial 
triad. We choose the orientation so that the magnetic field is aligned with e\, 
thus ensuring B2 = S3 = while B\ — h. From the Maxwell equations for B a 
this yields the constraints 

n 2 = -a 13 , n 3 = ct 12 (22) 

The fi 1 is still left undetermined, but to simplify the equations we will choose a 
frame so that fi 1 = 0. This choice of gauge makes the anisotropic-stress tensor 
diagonal (as in the general case of refs.jll]], [|l3)), with 

7r MV = D^h 2 (23) 

where D^ v = diag(0, — |, |, |). 

Introducing the parameters a±, defined by 

(oil, CT22, 033) = (-2o-+, cr+ + V3cr_, cr + - V3a-), 
allows the full set of equations to be written as a first-order dynamical system: 



d t h = — Bh-2a + h (24) 
3 

d tP = -jGp (25) 

d t e = -ie 2 -a 2 - l -[{z 1 -2) P + h 2 ] 

-|[4(3 7 - \)p 2 + (67 + A)ph 2 + 3/i 4 ] - 6 X U (26) 
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d t U 


4 

= 

O 


<V+ 


= -9(7+ 


d t cr- 


= -9(7- 




= -9(723 


dtcri2 


= -9(712 


d t <7l3 


= -9(713 


where \ — 


(k 2 A) _1 is c 
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+ ± a+ e 



'" + 3 h 



/ 2 

K 2 - ^ 



V3 

2(7l2CTl3 



making the redefinitions: nh — > h and up 
further symmetry in these equations. By a change of variables 



3ja +P (27) 

(3 7 -2)p+^ 2 ) (28) 

(29) 
(30) 

3^23 (31) 

2CT23, (32) 

ive set k — 1 implicitly by 
p. Interestingly, there is a 



0-12 4 

V3(7_ -4 

the equation for (f> turns into 



■1013 
2(723 



((7 B +Z(7 C )e 2 ^, 



9t</> = CFC 



(33) 
(34) 

(35) 



and can be completely eliminated from the system. Thus, this variable can be 
ignored. The remaining equations for the shear are 



d t <J+ = -9cr + - a\ + i/1 2 



1 



(37-2)p+-/i 2 



Sta^ = (-9 + 3(7+ + ob)oa 



d t <JB 
dt&c 



— 9(7 s — c 
(-9-2(7 B )(7 C 
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(36) 

(37) 
(38) 
(39) 



A further simplification is achieved by defining a new set of expansion nor- 
malized variables: 



3(7+ 

V3cTA,B,C 
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£+ = 

u = 
H = 
ft = 

We introduce a new time variable by 

dt _ 3 
d7 ~ 9 

and denote differentiation with respect to r by ', so the evolution equations 
become: 



92 

h 

9 
3p 

9 2 



(40) 
(41) 
(42) 
(43) 
(44) 

(45) 



9' 



(1 + 9)9 



(46) 
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u' = 

n' = 

W = 

E' + = 

5^ = 

X'b = 

E' = 



where 



E 
1 



Ti 2 

(2q — 2)17 — — (W + 5E+W 2 - 3 7 ft) 9 2 
18 

[2«-(37-2)]n 
(g-l-2E+)W 

(g - 2)E+ - 3Ei + 3W 2 (l - ~ [(37 - 2)fi + 4W 2 ] 6 : 

g - 2 + 3£+ + v/3£ B ) E A 
(g - 2)E B + 2y/SI% - ViY? A 
g-2-2\/3E B N ) E c 



(47) 

(48) 
(49) 

(50) 

(51) 
(52) 
(53) 



l + E + i(3 7 -4)ft 



+X 
Ei 



(37 

- E' 



n 2 



f ^ 7 ^ 2 + V 
4 4 

E 2 



X 



Q2 
6 



n 



-w 2 



e 2 



ft 4 



9 2 + 18x?7 



(54) 
(55) 
(56) 



and the last of these is the constraint equation. We note that in the limit of 
general relativity (x = 0), the variables U and decouple and f2 can then 
be determined from the constraint equation, leaving us with 5 fundamental 
variables: (H, £+, £a,b,c)- 



2.1 Flat Friedmann universe 

Let us first look at one of the simplest cases, the flat isotropic Friedmann brane- 
world. We set Ti = E + = Y,a,b,C = 0, and the solution can now be derived 
exactly. The result is a two parameter family (Qq,K) of exact solutions given 
by: 



fi( T ) = (|e^e~ 37T + 1 + Ke^~^ ' 
18 X U(t) = Ke^- i)T n(T) 

e(r) = e oe -^ T n(r)-^ 



(57) 

(58) 
(59) 



We note that if 7 < | the Friedmann brane-universe will be U dominated 
near the initial singularity. Otherwise, the quadratic matter-term will dominate 
initially in eqn. (|lj). This is quite different from the general relativistic case. 
In the general relativity case (\ = 0), the f2 term will dominate initially in 
an isotropic universe, where il = 1 exactly. The metric for a flat Friedmann 
brane-world can now be written as 



ds 2 



Q(ry 



dr 2 + e 2T (dx 2 + dy 2 +dz 2 ) 



(60) 



By integrating eqn. (Eq) we can restore the comoving proper time coordinate t. 
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2.2 Non-magnetic Bianchi type I solutions 

Now we examine the invariant solution subspace TL = ^a.c — of Bianchi 
type I universes containing perfect fluid and magnetic held. If we introduce an 
angular variable 4> and parameters A, K and Oo the solutions can be written in 
parametric form: 



18 X U(t 
6(r 
S + (r 



= Aei(T- 2 ) T fi(r)^cos0 
= Ae^t-^fliT^sincj). 



A 2 e 3( 7 -2)r 



(61) 

(62) 
(63) 
(64) 
(65) 



These solutions have been found earlier ||. Let us look at the initial sin- 
gularity in this case, which is approached as r — ► — oo. The behaviour near 
the singularity in this case depends strongly on the parameter 7 Q . If the per- 
fect fluid has 7 < 1 then the isotropic perfect fluid stresses are dominated by 
the shear anisotropy at early times and a Kasner-like singularity results (in the 
general relativistic case this occurs for 7 < 2 ), as r — > —00: 



n « A- 2 e 3 ^ T 
£+ ps cos <j> 
Y.g i=s sin <f>. 



(66) 
(67) 
(68) 



In this case, the behaviour is that of a vacuum Kasner universe, and the metric 
near the initial singularity can be approximated by the Kasner metric, 



ds 1 



-dt 2 



t^dx 2 +t 2 ^dy 2 +t 2p3 dz 2 



(69) 



where J^Pi — J2Pi = 1- ^ n the case of 7 = 1, the solutions can be written 
similar to the Kasner metric, but now the exponents fulfil the equations ^2pi — 
cosip > 0, YlPi — cos 2 -;/'- The solution space is thus as a half-sphere and the 
situation is analogous to that of general relativity when 7 = 2. 

If 7 > 1, we get a matter-dominated behaviour as r — > —00 (this type of 
behaviour cannot occur in general relativity because it would require 7 > 2) 
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xeg" 

VGA 

Vx&o 

V6A 



= 37T 



e (37-l)r COS( 



e (37-l)r sin< 



(70) 
(71) 

(72) 



In this case we see that the behaviour is quasi-isotropic. The shear is going 
exponentially towards zero in terms of r. The metric in this case can be ap- 
proximated by 



ds 2 



-dt 2 +t^ (e-^ t)coaW dx 2 +e^ (t)cos ( 0+ t)d 2/ 2 + e ^(*) cos (*-f)dz 2 ) (73) 



where 



This metric differs slightly from the Friedmann metric, since the quadratic mat- 
ter term dominates in the initial singularity. In some sense it asymptotes to a 
model with zero brane tension. It is a non-general-relativistic model which was 
first investigated by Binetruy, Deffayet and Langlois fl3l"| . 



3 The singularity in magnetic Bianchi type I 
brane-worlds 

We will now investigate a more general solution of the magnetic Bianchi type 
I brane-world. We have just seen that the quadratic density terms affect the 
singularity significantly in the non-magnetic case whenever 7 > 1. Therefore 
we expect that the anisotropic tracefree stresses contributed by the magnetic 
field will have a significant new anisotropising effect on the form of the solution 
near the singularity in many cases where 7 > 4/3. It is useful to compare 
with the general relativistic solutions with a magnetic field. If we set x = U = 
^A.c — SI = we get the so-called Rosen solutions. These are 2-parameter 
pure magnetic solutions, and an interesting point about these solutions is that 
the near the initial singularity the behaviour is that of a Kasner vacuum. The 
solutions can is mapped onto the straight line Sb = C(E + — 2) where we have 
the bounds El + E| < 1, 3C 2 < 1. 



3.1 Rosen branes 

Let us try to find similar solutions for the brane-world scenario. First we look 
at the invariant space Q = = E^.p. We make the ansatz 



(75) 



where p is a positive constant. Clearly TL oc e pT , thus TL is exponentially de- 
creasing towards the initial singularity r — > —00. Let us assume that we are in 
the vicinity of the Kasner circle, i.e. E + and E^ are approximately constants 
and q « 2. For this to hold TL 2 ® must also be exponentially decreasing towards 
the initial singularity. Since 0' = — (1 + q)Q, we get the inequality: 



q - 3 - 4E+ > => - E+ > - 



(76) 



In addition, we get El 
solutions of the form: 



T, 2 



= 1 and near the initial singularity there are 



e 
n 

E+ 
E s 



Hoe (l-2cos0)T 
COS(j) 

sin (b. 



(77) 
(78) 
(79) 
(80) 



9 



where the constant 4> has to satisfy 



cos0 < --. (81) 



The function U is 



Ue 2T + 2^. (§±§2S*V e -2(4co S ^+l)T 2C0S^^-1 

36 



(V 



2 cos d> = — 1 



3.2 Shear-dominated singularity with perfect fluid 

Let us again assume that the shear dominates near the initial singularity, but 
now we insert a perfect fluid into the equations, so we will only assume that 
Y,a.c = 0. We will also assume that we are approximately on the Kasner circle, 
i.e. E = 1. Thus we get q = 2 and E' = 0. The same reasoning as before yields 
the bound 

- E+ > 1 (83) 
In addition, from the Tt 2 £lQ 2 term, we get the bound 

-E+>~(3 7 -2). (84) 



We note that for 7 < 1 the bound ( J8 3| ) is sufficient for bound (|8J|) to hold. In the 
case 7 > 1 the bound (|84|) is sufficient for bound ( |S3"| ) to hold. The behaviour 
for O, 7i, E + and E# is the same to lowest order as in the previous example, 



except that E + has to satisfy both (83) and (|34j). The evolution for fl is to the 
lowest order 

Q = n e 3{2 -~< )T (85) 



3.3 Fluid-dominated singularity 

As we saw in the simple case with no magnetic field present, perfect-fluid matter 
can dominate the singularity if 7 > 1. This is also the case with a magnetic 
field. In that case we have the term <|£1 2 2 w 1. Let us assume that this is the 
case here. For this to happen we clearly require 

q = 3 7 -i. (86) 

On the other hand, the shear has to be exponentially decreasing in r as we go 
towards the initial singularity, so 

q-2>0 =S> 7 > 1 (87) 

But from the 7i 2 f20 2 term, we get an even stronger bound: 

7 > \ (88) 

The physical explanation for this bound in order that the solutions to be stable 
in the past is evident: if the perfect fluid dominates the past singularity, it 
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inevitably drives the universe towards isotropy. Near isotropy the magnetic 
field behaves as a ^ EM — | perfect fluid, thus if 7 < |, the magnetic field will 
become dominant in the past. The fluid-dominated universe will therefore be 
unstable in the past for 7 < |. If the inequality, 7 > |, is fulfilled then the 
evolution near the initial singularity is approximated by 

6 = 9 e- 37T , (89) 

n = -=|U 37T , (90) 

H = W e (37 - 2)T , (91) 
S = S e 6(7 ~ 1)T . (92) 

3.4 ^/-dominated and 7i-dominated singularity 

Let us look at the other terms in the constraint equation and check whether 
the dark energy, U, or the magnetic field, H, can be dominant near the initial 
singularity. As we have already seen, the U term dominates in the Friedmann 
case if 7 < | . This is the only situation where the U term is allowed to dominate. 
For the y/xH 2< d = H term, the constraint equation allows this term to grow to 
arbitrary large values. Let us consider the case where H diverges. The equation 
for H is 

H 1 = (q - 3 - 4E+)H. (93) 

For H to diverge at the initial singularity, (q — 3 — 4E + ) < 0. This is not possible 
if H is large. Thus H cannot diverge in the initial singularity. 
We note however that there are solutions of the form 

H 2 ex — (94) 

K - q T 

£+ ex — (95) 

K -q r 

but these do not describe a physical singularity. These are merely a result of 
the choice of time coordinate. These solutions yield 

9 oc K - q r, (96) 

and are only turning points in the evolution of the brane, where the brane 
recollapses. The infinity in H in the phase space does not represent a physical 
singularity in the brane evolution. 

There are also solutions where H is initially a constant. One of these has a 
particular value of q: 



1= ^(89 + V3345) (97) 

The solutions in this case are 

£+ = ^(23 + V3345) (98) 
08 

9 = 6 e- 4 ( 1+s +) r (99) 
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1(2 -37) - 1 - 



4 



Figure 1: The stability of the Kasner circle. The positive signs mean stability 
in the past for the different eigenvalues. 

S A b,c = £ 0: AB,ce (1+4S+)T (100) 
H = W e 2(1+E+)T (101) 
n = ft e (4+8E+ ~ 37)r (102) 

This solution corresponds to a stable fixed point, but is not an exact solution 
because H = implies H = 0. The fixed point correspond to a non general- 
relativistic model, with the brane tension zero. Writing out the asymptotic form 
of the metric in the past 

ds 2 = -dt 2 + i-M(v^345-43) da .2 + t ^ dy 2 + (1Q3) 

The magnetic field diverges like 

h * = J^Lt- 1 (104) 

1 + 9 

where q is given as above, and H = ^v/2305 + 44V3345. Thus in this particular 
model, there are no free parameters prescribing the metric or the magnetic field. 
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3.5 (In) Stability of the Kasner circle 



The (in) stability of the Kasner circle proved to be an important issue in the 
general relativistic case M] . We have already noticed that the Kasner solutions 
are unstable for E + > — jj~But we can show that all points on the Kasner circle 
are saddle points. In our derivations we have put £U,c = 0. Linearising the 
equations for Ea,c around the Kasner solutions yields 

Y,' A = (3£+ + \/3£b)£a, (105) 
T,' c = -2V3£bI] c . (106) 

Thus from this alone we can say that each point on the Kasner circle is a 
saddle point for past evolution. Therefore the generic solution does not "start" 
at the Kasner solution. The stability analysis of the Kasner circle is illustrated 
in Figure [l]. The two horizontal lines arise from the TL 2 Q 2 £l and the 7i 4 2 
terms. 



4 Equilibrium points and stability 

Due to the complexity and the non-linearity of the equations the we will intro- 
duce the following new variables: 

h = Vxn 2 e (107) 
lo = y/xne (108) 

V = 18 X U (109) 

r = n 2 (no) 

but we will still keep the variable f2. The cost is that we get another constraint: 

HVl = uoR (111) 

The variable V can be in principle be obtained from the constraint equation. 

The set of equations now becomes simpler for stability analysis, but we see that 

part of the boundary of M = {(H, R ) ui,n > 'E +) Y<a,b,c)\HVL = ujR and H,R,Lo,fl>0} 

has to be excluded. The boundary is important, however, because many of the 

equilibrium points lie on the boundary. We will therefore include all of the 

boundary in the analysis, but bear in mind that part of the boundary is un- 

physical. 



We must also note that the hypersurface given by (111) has a conical singu- 
larity at the origin. The tangent space at the origin therefore has a degeneracy, 
(in this case is a U (1) x U (1) degeneracy). This must be taken into account when 
performing a perturbation analysis around these singular points, but it is advan- 
tageous that the singular space is an invariant subspace of the dynamical system. 
Hence, through every point p 6 M there is a unique maximally-extended orbit, 
passing through p. 

We consider the set M. parametrised by (H, R, lo, Q, E+, £a,b,c)- The general- 
relativistic case is defined by the invariant subset H = lo = V = 0, but if 
H = lo = and V > 0, then V can be interpreted as a radiation fluid in the 
general relativity case (note that in general V can be negative). Let us exam- 
ine the interesting equilibrium points. All of those with lo = H = V — are 
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general-relativistic equilibrium points and were previously found in rcf . |T6| . We 
check whether their stability properties change after the inclusion of the "non- 
general relativistic" variables u>, H and V. We will also carry out a stability 
analysis for the new equilibrium points. The results are as follows: 

1. P(I), flat Friedmann universe 

R = Y, + = Y, ABC = H = u = V = 0. 
g =|(3 7 -2),'n = l ) 0< 7 <2. 

This is a future attractor in the general relativistic models when < 7 < | . 
It is also future stable in the H and u> direction for < 7 < | and thus 
unstable to the past. 

2. CR(I) flat Friedmann universe with radiation fluid. 
R = S + = Y.a.b,c = H = uj = 

q= 1, fi + y = 1, 7= f, Q> 0. 

This is a line bifurcation which is a future attractor. This set of points 
connects the P(I) and PV(I) equilibrium points. 

3. PV(I), flat Friedmann universe with V fluid. 
R = £+ = Y, a ,b,c = H = u =O = 

q= 1, V = 1, < 7 < 2. 

This equilibrium point is similar to that of a radiation-dominated universe, 
and is stable in the future for 7 > |. 

4. PQiOO 

i? = £+ = Sa.s.c = # = ^ = ^ = 
g = 37-l, uj'='VE, 0<7<2 

This model is an attractor in the past for 7 > |. In the shear variables 
the model is a past attractor for 7 > 1, but the magnetic variable H is 
not stable in the past for 7 < |. Thus, for the generic solutions, this is a 
saddle point in the past for 7 < |. 

5. £Q 2 (J) 

H = S + = S A!B> c = J ff = n = 

This is a one-parameter bifurcation. Each point on it is a saddle point. 
It serves as a one-parameter set of stable points to the future for £1 = 0. 
However, this is not physically possible in the original model, since f2 = 
implies u> = 0. 

6. flat Friedmann universe with vacuum fluid 
R = £+ = £ab,c = H = V = 

q = -1, 7 = 0, £ + = 1 

This is a one-parameter bifurcation that is an attractor for future solu- 
tions. 

7. PMi(J) 

£4.b,c = ^ = P = l/ = 
fl=|(2- 7 )(37-4), S+=i(3 7 -4), 
flf=|(37-2),n = |(4- 7 ),|<7<2 

In the general relativistic case this is an attractor in the future for | < 
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7 < I , but becomes a saddle point in the brane- world model because it is 
unstable in the V-direction for 7 > |. 

8. PM 2 (I) 

Z c = uj = H = V = 

R=\{2- 7) (97 - 14), £+ = |(3 7 - 4), £4 = v /|(2- 7 )(5 7 -8), 

S B = -iV3(5 7 -8), 

9=1(37-2), fi = |(2-7), f<7< 2 

In the general relativistic case this is an attractor in the future for | < 
7 < I , but becomes a saddle point in the brane- world model because it is 
unstable in the Indirection for 7 > |. 

9. PM 3 (I) 

lu = H = Q = V = 

r> 1 v 1 V / 7 v 1 V — 1 

«- 3' ^+ - 4> ^ - y 24, - - 

?=|, 0< 7 <2 

In the general relativistic case this is an attractor in the future for | < 
7 < 2, but becomes a saddle point in the brane- world model because it is 
unstable in the ^-direction. 

10. JC, Kasner vacuum 

R = £ A ,c = Q = u> = H = V = 
£^_ + £| = 1, q = 2, < 7 < 2 

As we have seen, the Kasner circle is a one-parameter set of points, all are 
saddle points for the generic brane- world solutions. 

11. SIC, Kasner-dust half-sphere 
R = £ A ,c = ft = H = V = 
£2 + £ = l, g = 2, 7 = 1 

This is a two-parameter bifurcation set. The whole set is a set of saddle 
points for the system of equations. 

12. CM (I) 

H = lu = V = 



R — V2 1 1 V _ 1 V — 1 /q v2 - 1 v2 1 l^v2^7 

' 24' Zj + 4' I^V 13 ' — 2^.4 48' 24 A ^ 24' 

« - 2 O — T _ *v — ^ 

1 — 2' " — 8 A' 7— 3 



In the general relativistic case this is a one-parameter line-bifurcation 
that attracts future solutions, but in the brane-world model the points 
are unstable in the V direction. 

13. Pffi(I) 

R = lu = Q = £a,b,c = 

s+ = 23+^3345 ) H 2 = i s+(1 + 4S+)) y = ^(13725 + 107^3345), 
q = ^(89 + V3345), < 7 < 2 

This is a new and unusual fixed point, and is an attractor to the past. 

14. PH 2 (I) 

R = Lo = n = £4,5,0 = 

£+ = 23 '^ , H 2 = i£+(l + 4£+), V = ^i( 13725 - 107V3345), 
q = ^(89 - V3345), < 7 < 2 

This is similar to the fixed point PM\(I) but it is a saddle point. 
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15. 



R = f2 = T^a,b.c — 

S + = i(3 7 -4) )g = 3 7 -l 

H 2 



j_ 

2a 



Vb 2 



7 e [i(lll-x/3345),l 
where a 



4oc), a; 
U 



|(7-l)(3 7 -4)-4g 2 
H(3 7 -2) 



4 

45 7 -22, 6 = 3 (^- 7 3 - 114 7 2 + 146 7 - 56)and c = §(7 - 1) 2 (3 7 
The local behaviour around this fixed-point turns out to be quite complex, 
but a careful analysis shows that this fixed point is saddle point for all 
allowed values of 7. For 7 — ^(111 — V3345), 1 and 4/3, H = and will 
coincide with one of the points of PH2{I),SK. and PQ\(I) respectively. 

16. PN 2 {I) 

R = ft = Sa.b.c = 



[§,2] 



E + = ±(3 7 -4),< z = 3 7 - 
ff 2 



4' 
J_ 

2« 



-6 + \lb 2 - 4ac) , w = 



k7-l)(37-4)-4g 2 
ff(3 7 -2) 



45 ' 3 



> 



where a, b and c are given in PNi(I). This point is a saddle point for all 
7- 

We can note one particular thing, the properties of the future attractor points 
are altered from the general relativistic case. In the brane- world model the dark- 
radiation term will make the general-relativistic future attractors unstable for 
7 > |. Thus, both the local and global properties of the phase space are altered 
on passing from the general relativistic to the brane-world model. For example, 
we have brane-world solutions that fly off to infinity in finite time while in 
the general relativistic case the phase space is compact. Also, the brane-world 
solutions behave differently in the past. The general relativistic case had no 
attractors in the past for any 7. But as we can see, in the brane-world model, 
there are two past attractors for 7 > | (PQi(I) and PHi(I)) and one for 7 < | 
(l'Ih(l))- 

The importance of these fixed points can be expressed in the case 7 > ^ by 
the following theorem: 



Theore m Let p be a point in the configuration space M. satisfying the con- 



straint (111), and let s p (t) be the maximally extended solution satisfying the 
equations of motion such that s p (0) = p. Then for 7 > | there exists a compact 
region C C M and a t such that for every t < tq then s p (t) € C. 



V 2 
A 



Proof: Consider the compact region T C M. given by E = E 2 
T? c < E* > 1, Q < fl* > 0, uj < uj* > V6, R < R* > 0, H < H* > 
assume that 7 > |. For E > 1 we have E > |E + | and 

E' = 2(q - 2)E + 6i?E + - [(3 7 

E 



= 2(,7-2)E + 6i?E + 

> 2(E — 1)E + 6i?E_| 

> 2(E-l-6i?)E 



2)ujH + AH 2 } E_ 



-H 2 - i(3 7 - A)ojH + -uj 2 
2 2 V ' ' 3 



E| + 
I and 



(112) 
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/C, Sj 



PHiiX) 




Figure 2: A flow-diagram for the non-general-relativistic H — £+ plane where 
R = u> = Q = T,a.b,c = 0. Five equilibrium points can be seen; four are marked 
on the figure. The last one PH 2 (I), is a saddle point and can be seen to the 
left and below the centre of the figure. 
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Figure 3: A diagram of the web of equilibrium points. Horizontal lines are 
bifurcation sets; all are line bifurcations except for SIC, which is a half-sphere. 
Lines in red are future attractors, while lines in green are past attractors. 
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because the polynomial in lo and H inside the square brackets is positive definite 
for 7 E [§, 2]. Similarly, for H > H* we have 



#'> 12 ^H"*- y-J (113) 

and for u) > lo* we have 

u/>(3 7 -2)(^-l)u;* (114) 

These equations ensures that the variables H and lo are finite in the past. We 
also have for H ^ 

i)'-^" (#)*'(#) <115) 

Integrating yields the bound j|Q < -g^ye 27 " for r < 0. Hence, in the past R 
can be driven arbitrarily close to zero as r — > — oo since -ff is bounded. Thus 
for S > £*, £' > for some time in the past. Hence, £ is bounded in the past 
too. Similarly, for lo ^ 0, ft can also be driven arbitrarily close to zero. Thus in 
this case, the region T is sufficient and we can set C = T . 

If H = (which is an invariant solution subspace) the equation for V is 

V' = (2q-2)V (116) 

which after a rearranging yields 

(ln|F|)' = 2g-2>0 (117) 

Thus V is bounded in the past, and according to the constraint equation all of 
the other variables have to be as well. 

If w = but H > 0, then fT > (37- 4) (ft* - l)fi* > for ft > ft* > 1 which 
makes ft bounded in the past. 

The theorem now follows from this analysis. I 



This theorem shows that the solution cannot fly off to infinity as r — > —00. 
With some amendments, we believe that the behaviour at the initial singularity 
is "well-behaved" for any 7, even though we have not been able to be prove a 
similar theorem in that case. 

In Figure |^ we have drawn a flow-diagram of the dynamical system. The 
section shown is a part of the non-physical boundary, and the vector field is 
projected onto this plane. In Figure |^ the web of equilibrium points are illus- 
trated. We see how the points are connected through bifurcation sets, and it is 
clear how some branch off from other equilibrium points. 

From this analysis we believe that it is very likely that the pure magnetic 
Bianchi type I brane- worlds do not exhibit chaotic behaviour in the past. 



5 Discussion 

In this paper we have investigated pure magnetic Bianchi type I brane- worlds. 
We found asymptotic forms for the solutions near the initial singularity for 
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different configurations of matter, magnetic field and shear anisotropy.. We also 
checked whether there are any equilibrium points that could serve as a past 
attractor for the generic solutions. For 7 > | we found two such attractors, 
while in the case 7 < | there was one. This may indicate that there can 
be no chaotic behaviour in the generic solutions for a magnetic Bianchi type I 
brane- world. To the past, the solution can approach the general relativistic case 
but will only oscillate a finite number of times. The possible initial states of 
these two equilibrium points, correspond to two different non general-relativistic 
models, one of them appears to be previously unknown. 

When we investigated the various equilibrium points, we noted that all of 
the future stable general-relativistic points were unstable in the brane-world 
scenario for 7 > |. Another major difference between the two cases is that 
the configuration space is non-compact for the brane-world solutions. We found 
solutions that fled to infinity in a finite coordinate time, and found that the non- 
local energy 11 can become arbitrary large (and even negative) and so presents 
a major obstacle to this sort of analysis. 

The models examined in this paper are based upon the simplest anisotropic 
universes. It remains to be seen what new features might emerge from a study 
of the general case. However, we know that in the magnetic case (in contrast 
to the non-magnetic situation) many of the features of the general anisotropic 
universe are present in the simple type I model. 

Note: After this paper was completed and submitted, a similar and indepen- 
dent work by Alan Coley appeared |52], |33| addressing some of the same issues 
as we have done in this paper. 



Acknowledgments 

SH was funded by The Research Council of Norway. 

References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 



L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) 
L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) 
T. Shiromizu, K.I. Maeda, and M. Sasaki, Phys. Rev. D 62, 024012 (2000) 



J.D. Barrow and R. Maartens, ;r-qc/0 108073 



A. Campos and C. F. Sopuerta, Phys. Rev. D 63, 104012 (2001) 
A. V. Toporensky, Class. Quant. Grav. 18, 2311 (2001) 



M. G. Santos, F. Vernizzi and P. G. Ferreira, hcp-ph/0103112. 



A. Campos and C. F. Sopuerta, istro-ph/0 102355 



C. Chen, T. Harko and M. K. Mak, hcp-th/0103240 



R. Maartens, V. Sahni and T.D. Saini, Phys. Rev. D 63, 063509 (2001) 



20 



[11] J. D. Barrow, Phys. Rev. D 55, 7451 (1997) 

[12] J.D. Barrow, Phys. Rev. D51, 3113 (1995) 

[13] J.D. Barrow and R. Maartens, Phys. Rev. D 59, 043502, (1998) 

[14] V. A. Ruban, Sov. Astron. 26 632 (1982) 

[15] S. Hervik, Class. Quant. Grav. 17, 2765 (2000) 

[16] V. G. LeBlanc, Class. Quantum Grav. 14, 2281 (1997) 

[17] B. L. Spokoiny, Phys. Lett. A 81, 493 (1981) 

[18] V. A. Belinsky, I. M. Khalatnikov and E. M. Lifshitz, Adv. Phys. 19, 525 
(1970) 

[19] J.D. Barrow, Phys. Rev. Lett. 46, 963 (1981) 

[20] J.D. Barrow, Phys. Reports 85, 1 (1982) 

[21] D. Chernoff and J.D. Barrow, Phys. Rev. Lett. 50, 134 (1983) 

[22] J. Demaret, M. Henneaux, and P. Spindel, Phys. Lett. B 164, 27 (1985) 

[23] T. Damour, M. Henneaux, B. Julia and H. Nicolai, Phys. Lett. B 509, 323 
(2001) 

[24] T. Damour and M. Henneaux, Phys.Rev.Lett. 86, 4749-4752 (2001) 

[25] J.D. Barrow and K. Kunze, Phys. Rev. D 55, 623 (1997) 

[26] C. Tsagas and J.D. Barrow, Class. Quantum Gravity 15, 3523 (1998) 

[27] C. Tsagas and J.D. Barrow, Class. Quantum Gravity 14, 2539 (1997) 

[28] L. P. Hughston and K. J. Jacobs, Astrophys. J. 160, 147 (1970) 



[29] R. Maartens, fer-qc/0101059 



[30] R. Maartens, Phys. Rev. D 62, 084023 (2000) 

[31] P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B565, 269 (2000) 
[32] A. Coley, |hep-th/0110049| 



[33] A. Coley, bep-th/0110117 



21 



